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We report a study of the decay B° — >• DK + n~ followed by D K~n + , where D indicates D° 
or D° . We reconstruct the DK + ty~ state in a phase space corresponding to D_K'*(892) . The CP- 
violating angle 03 affects its decay rate via the interference between b — >• u and b — > c transitions. 
The result is obtained from a 711 fb _1 data sample that contains 772 x 10 6 BB pairs collected at the 
T(4S) resonance with the Belle detector at the KEKB asymmetric-energy e + e~ collider. We measure 
the ratio TZ DK , a = r(B° ->■ [K~-k+] d K +tx~)/T{B° -> [K+-r-] d K +tt~) to be (4.1±^g±?-.|) x HT 2 , 
and set an upper limit of 1Z DK *o < 0.16 at the 95% confidence level. 



PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Hh, 14.40. Nd 



Determination of the parameters of the standard model 
is important as a consistency check and as a way to 
search for new physics. In the standard model, the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix [lj V con- 
sists of four independent weak interaction parameters 
for the quark sector; the three CP- violating phases 0i, 

02 and 03 are defined as the angles of one particu- 
lar CKM unitarity triangle with the latter defined as 

03 = a,Tg(-V u dV u b* /V cd V cb *). This phase is less accu- 
rately determined than the other two 0. In the usual 
quark-phase convention where large complex phases ap- 
pear only in V u b and Vtd 0, the measurement of 03 is 
equivalent to the extraction of the phase of V u b relative to 
the phases of other CKM matrix elements. To date, the 
03 measurement has been advanced mainly by exploi ting 
charged B meson decays into D^K + final states [44l2| 
wherein the CP sensitivity is due to the interference be- 
tween the two amplitudes of Z)(*)° and decays into 
a common final state. 



In this paper, we consider the neutral meson decay 
B° — > DK*° as an alternative process for measuring the 
angle 03 . As shown by the Feynman diagrams in Fig. [1] a 
weak decay of the B meson is tagged by the K*° decaying 
into K + tt~ [l3| . We measure the ratio 'R- DK *o (Til. HEf 



B° 



D" 



B° 



FIG. 1: Diagrams for the B° -> D°K*° and B° -> D°K*° 
decays. The 03 dependence in the b — > u transition is ex- 
tracted from the interference of the two decay paths, which 
occurs when the D and D° mesons decay to the same final 
state. 



defined as 



TZ DK *o = 



T(B° -> [K-tt+} d K+tt- 



T(B° -> [K+-k-] d K+-k-) 
r s + r D + 2kr s r D cos (6 S 



S D )cos 03,(1) 



where r D = \A(D° -> K +ir-)/A{D° -> R-Tr+)\ is the 
ratio for D decay amplitudes and 5d is the strong phase 
difference of the two D decays appearing in this ratio. 
Both r_o and Sd have been obtained experimentally [lij . 
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The parameters rs, 6s and k are defined as 

= r(B°^D°K+K-) = JdpA^ u (p) 
s F(B° -> D°K+it-) fd P A 2 b ^ c (p)' [) 

where Ab_>. c (p) and At,_). u (p) are the magnitudes of the 
amplitudes for the b — > c and b — > u transitions, respec- 
tively, and S(p) is the relative strong phase. The variable 
p indicates the position in the DK + ir~ Dalitz plot. In 
this analysis, we calculate the integrals over a phase space 
of the state DK*(892)°. In the case of a two-body B de- 
cay, rs becomes the ratio of the amplitudes for b — > u and 
b — > c and k becomes 1. The value of rs is expected to 
be around 0.4, which is obtained from |KbV c * |/|V c bV^*J 
and depends on strong interaction effects. According to 
a simulation study using a Dalitz model based on recent 
measurements [17], the value of k is around 0.95 in the 
phase space of interest here. One observable IZdk* is 
not enough to extract the four unknowns fa, rs, k, and 
5s- However, the measurements for other D decays such 
as D — >• K + K~ and Ksir provide additional informa- 
tion needed to extract fa, where the observable TZrjK* 
should be defined in the same phase space of the B° decay 
between different D decays so that the same parameters 
rs, k, and 6s can be used. The decay in the numerator 
of Eq. ([I]) is the signal mode, referred to as the "sup- 
pressed mode," while the decay in the denominator is 
the calibration mode referred to as the "favored mode." 

This result is based on a data sample that contains 
772 xlO 6 BB pairs, collected with the Belle detector 
at the KEKB asymmetric-energy e + e~ (3.5 on 8 GeV) 
collider [l8| operating at the T(45) resonance. The 
Belle detector is a large-solid-angle magnetic spectrom- 
eter that consists of a silicon vertex detector, a 50-layer 
central drift chamber (CDC), an array of aerogel thresh- 
old Cherenkov counters (ACC), a barrel-like arrangement 
of time-of- flight scintillation counters (TOF), and an elec- 
tromagnetic calorimeter comprised of CsI(Tl) crystals lo- 
cated inside a superconducting solenoid coil that provides 
a 1.5 T magnetic field. An iron flux- return located out- 
side of the coil is instrumented to detect K\ mesons and 
to identify muons. The detector is described in detail 
elsewhere [19j . 

Charged kaon and pion candidates are identified us- 
ing ionization loss in the CDC and information from the 
ACC and the TOF. The efficiency is 85-95% and the 
probability of misidentification is 10-20%. We recon- 
struct D mesons from pairs of oppositely-charged kaon 
and pion candidates. We require that the invariant mass 
is within ±15 McV/c 2 (±3cr) of the nominal D° mass. 
K*° candidates are reconstructed from K + ir~ pairs. We 
require that the invariant mass is within ±50 MeV/c 2 of 
the nominal K*° mass. We combine D and K*° candi- 



dates to form B° mesons. Candidate events are identi- 
fied by the energy difference AE = ^ ■ E{ — Eh and the 
beam-constrained mass M^ c = \J E^ — J2iPi\ 2 > where 
i?b is the beam energy and pi and Ei are the momenta 
and energies, respectively, of the B° meson decay prod- 
ucts in the e + e~ center-of-mass (CM) frame. We select 
events with 5.271 GeV/c 2 < M bc < 5.287 GeV/c 2 and 
-0.1 GeV < AE < 0.3 GeV. In the rare case where 
there are multiple candidates in an event, the candidate 
with Mbc closest to its nominal value is chosen. 

Among other B decays, the most serious back- 
ground for the suppressed mode comes from B° — > 
[K*°K + ] D +ir~ . This decay produces the same final state 
as the B° —> DK* a signal, and the product branching 
fraction is about 10 times higher than that expected for 
the signal. To suppress this background, we exclude can- 
didates for which the invariant mass of the K ~ir + K + 
system is within ±18 MeV/c 2 (±3cr) of the nominal D + 
mass. The relative loss in the signal efficiency is 0.5%. 

Large combinatorial background of true D° and ran- 
dom K + and tt~ combinations from the e + e~ —¥ cc pro- 
cess and other BB decays is reduced if the D° is a decay 
product of D* + — > D Q ir + by using the mass difference 
AM between the [K^tt + ]dtt + and [K~it + ]d systems, 
where a tt + candidate is added to the latter to form the 
former. If AM > 0.15 GeV/c 2 for any additional 7r+ 
candidate not used in the B candidate reconstruction, 
the event is retained. This requirement removes 24% of 
cc background and 14% of BB background according to 
Monte Carlo (MC) simulation. The relative loss in signal 
efficiency is 5.0%. 

To discriminate the large combinatorial background 
dominated by the two-jet-like e + e~ — > qq continuum pro- 
cess, where q indicates u, d, s or c, a multivariate anal- 
ysis is performed using the following nine variables. 1) 
A variable obtained from the Fisher discriminants based 
on modified Fox- Wolfram moments [2(| where the coef- 
ficients of the Fisher discriminants are optimized using 
the signal and qq MC samples. This variable exploits the 
event topology, which is spherical and jet-like for BB and 
qq events, respectively. 2) The angle in the CM frame be- 
tween the thrust axes of the B decay and the detected 
remainders. For the latter, we assign the pion mass to all 
the charged particles and use photons with energy above 
0.1 GeV. 3) The signed difference of the vertices between 
the B candidate and the remaining charged tracks. For 
the signal event, the absolute value tends to be larger 
because of the longer lifetime of the B meson. 4) The 
angle between the K candidate from the D decay and 
the B candidate in the rest frame of the D candidate. 
Its distribution is flat for signal events but peaked near 
the extreme values for qq background. 5) The expected 
flavor dilution factor described in Ref. [2lJ. It ranges 
from zero for no flavor information to unity for unam- 
biguous flavor assignment. B candidates tend to have a 
larger flavor dilution factor than qq background. 6) The 
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angle 9 between the B meson momentum direction and 
the beam axis in the CM frame. The B decays follow a 
1 — cos 2 9 distribution, while the qq background is nearly 
flat in cos 9. 7) The distance of closest approach between 
the trajectories of the K* and D candidates. The value is 
close to zero for the signal but tends to be larger for the 
cc background. 8) The difference between the sum of the 
particle charges in the D hemisphere and the sum in the 
opposite hemisphere, excluding those used in the recon- 
struction of the B meson. The average charge difference 
is for the signal events but ±4/3 for the cc events, de- 
pending on the flavor of the B candidate. 9) The angle 
between the D and T(45) directions in the rest frame of 
the B candidate. The cosine distribution is about flat for 
signal events but peaks toward +1 for cc events. 

To effectively combine these nine variables, we em- 
ploy the NeuroBayes neural network package [22]. The 
NeuroBayes output is denoted as Cnb with a range of 
[— 1, f]. For example, events at Cnb ~ 1 are signal- 
like and events at Cnb ~ — 1 are qq-like. The training 
for the neural network optimization is performed by us- 
ing the signal and the qq MC samples, each of which 
contains 1 00,000 events after the event-selection require- 
ments. For the latter sample, we loosen the requirement 
on M bc to 5.23 GeV/c 2 < M hc < 5.27 GeV/c 2 to obtain 
a larger number of events, since all the input parameters 
have little correlation with Mb c - 

The Cnb distribution peaks at | Cnb | ~ 1 and is there- 
fore difficult to represent with a simple analytic function. 
However, the transformed variable 
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(4) 



where Cnb, low = —0.6 and Cnb, high = 1.0, has a distri- 
bution that can be modelled by a Gaussian. The events 
with Cnb < —0.6 are rejected. The background rejection 
rate is 70.5%, while the signal loss is 3.9%. 

The number of signal events is obtained by a two- 
dimensional unbinned extended maximum likelihood fit 
to AE and C^ B . The fits are applied separately for fa- 
vored and suppressed modes. For both modes, we catego- 
rize five common contributions. These are the DK*° sig- 
nal, the D°p° background, the combinatorial BB back- 
ground, the qq background, and the backgrounds that 
have peaks in the signal region of AE and C^ B ("peak- 
ing background"). In the favored mode, we include two 
more components: D°K+ and D°tt+. The B° -> D°p° 
decay satisfies the selection criteria when a pion from 
the p° decay is misidentified as a kaon. This compo- 
nent also includes other decays that satisfy the selec- 
tion criteria when a pion in the final state is misiden- 
tified as a kaon. The peaking background for the sup- 
pressed mode consists of B° — > [K~ir + ir~] D - K + and 
B° — > [K + K~] D ait + ix~ while the peaking background 
for the favored mode consists of B° — > [K + ir~ir~] D - K + . 
For the B° —> D°K + and D Q ir + backgrounds, a pion 



candidate is added to reconstruct K*°, where the lat- 
ter satisfies the selection when the ir + is misidentified 
as K + . We prepare two-dimensional probability density 
functions (PDFs) for each component as a product of 
one-dimensional PDFs on AE and C^ B , since the corre- 
lation between AE and C^ B is found to be small. 

The AE PDFs for a favored mode are parameterized 
by a double Gaussian for signal, a double Gaussian for 
D°p°, an exponential function for BB background, a lin- 
ear function for qq background, a Crystal Ball function 
for D°K + , and a double bifurcated Gaussian for D°ir + . 
The means and widths of the double-Gaussian PDFs for 
the signal and D°p° components are fixed from MC sam- 
ples. The mean of the AE distribution for D°p° is higher 
than that for the signal by about 70 MeV due to misiden- 
tification of a pion as a kaon. The parameters of the ex- 
ponential and linear PDFs are allowed to float. The AE 
PDF for the peaking background is defined to be that 
of the signal, and the yield is fixed by the world-average 
value of the branching fraction ]23]. The mean values 
of AE for D°K + and D°ir + are higher than those for 
the signal due to one additional pion and a misidentifi- 
cation for the latter mode. The shape parameters of the 
AE PDFs for these components are determined from MC 
and their yields are fixed by the world-average value of 
the branching fraction. 

The C^ B PDF is a sum of two Gaussians for each com- 
ponent. The shapes for the signal and BB background 
are fixed from the MC samples of each decay model. The 
C NB PDF for D Q p Q is defined by the same function as 
that of BB background. The C NB PDF for the peaking 
background is described as a weighted sum of MC-based 
PDFs for all the constituents. The shape for the qq back- 
ground is fixed from the M^ c sideband data sample in 
data defined by 5.23 GeV/c 2 < M bc < 5.27 GeV/c 2 . The 
validity of this use of the M^ c sideband sample, which is 
reasonable since all inputs for C NB have little correla- 
tion with Mbc, is checked using MC samples. The C^, B 
PDF for D°K+ and Dir+ is the same as that of the BB 
background. 

The results of the fits for suppressed and favored modes 
are shown in Fig. [2] and presented in Table Q] We obtain 
the ratio IZ^^' to be 



n 



where N< 



DK *° ' AT I 

JVfav/Cfav 

= (^lt±li) x io- 2 , 

sup (fav) is the signal yield for the suppressed (fa- 
vored) mode, e sup (f av ) is the detection efficiency obtained 
from a MC study for the suppressed (favored) mode. 

We list the sources of systematic uncertainties in Ta- 
blefll] The uncertainties of the PDF shape parameters are 
estimated by varying the determined parameters of the 
PDFs by ±l(j. The uncertainties due to the C NB PDFs 
for D°p°, combinatorial BB, D°K + , and D°ir + are esti- 
mated by replacing their PDFs with the signal PDF. The 
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FIG. 2: The projections of the fits to data for the suppressed 
mode (upper) and the favored mode (lower): the AE pro- 
jection for 3 < C*nb < 10 (left) and the C^b projection 
for \AE\ < 0.03 GeV (right). The fitted data samples are 
shown by the dots with error bars and the total PDFs are 
shown by the solid blue curve. Individual components are 
shown by the dashed red (DK*° signal), the dash-dotted 
magenta (D°p°), the short dashed green (combinatorial BB 
background), the long dashed brown (qq background), the 
very long dashed black (peaking backgrounds), the dash- 
dot-dotted gray (D°K + ), and the dash-dot-dot-dotted aqua 

(D°TT + ). 



uncertainty due to the PDF shape for qq is the largest sys- 
tematic uncertainty. The uncertainty due to the yields 
of the peaking background is conservatively estimated 
by applying and 2 times the nominal expected yields. 
The systematic uncertainty associated with the peaking 
background is small because of its small expected yield. 
The uncertainty due to the yields of D°K + and D°tt + 
is estimated by taking into account the uncertainty of 
the efficiencies and the branching fractions. We check 
the fit bias by generating 10,000 pseudo-experiments for 
each of the suppressed and favored modes. We obtain 
an almost standard Gaussian distribution for the pull, 
and take the product of the mean of the pull and the 
error of the nominal fit. MC statistics and the uncertain- 
ties in the efficiencies of particle identification dominate 
the systematic uncertainty in detection efficiency. The 
uncertainties in the efficiencies of particle identifications 
are determined from the decay D* + — > D°ir + followed 
by D° — >■ K~ir + . The uncertainty due to the charmless 
B° — > K*°K + tt~ decay is obtained from the upper limit 
of its branching ratio [23| and the efficiency estimated 



TABLE I: Summary of the results. The errors for TV and 
Tt DK ,o are statistical only. 





Mode 


e{%) 


N 


1l DK *0 


B°- 


■> [K~ it ]d K*° 


21.0 ± 0.3 
20.9 ±0.3 


r, v +10.6 
'■'-9.5 


(4.1^; 6 )xl0- 2 



TABLE II: Summary of the systematic uncertainties for 
TZ DK ,o. 



Source 


Uncertainty [10 2 ] 


Signal PDFs 


+0.1 - 0.2 


D°p° PDFs 


+0.0 - 0.1 


Combinatorial BB PDFs 


+1.8 - 1.2 


Peaking background PDFs 


+0.1 - 0.1 


qq PDFs 


+2.2 - 1.4 


D°K+ PDFs 


+0.0 - 0.0 


D n+ PDFs 


+0.0 - 0.1 


Fit bias 


+0.4 - 0.0 


Efficiency 


+0.1 - 0.1 


Charmless decay 


+0.0 - 0.3 


Total 


+2.8 - 1.8 



by assuming a non-resonant distribution in phase space. 
The uncertainties due to the favored mode are estimated 
in a similar manner as for the suppressed mode and are 
found to be small. 

The distribution of the likelihood C is obtained by con- 
volving the likelihood in the (AE,C^ B ) two-dimensional 
fit and an asymmetric Gaussian whose widths are the 
negative and positive systematic errors. We set a 95% 
confidence level (C.L.) upper limit for IZdk* to be 
IZok* < 0.16. We obtain an upper limit of rs < 0.4, 
corresponding to 95% C.L. limit of Rjjk* , by conserva- 
tively assuming that rs is much larger than rn so that 
Rdk'° = r s- The uncertainties due to the signal yield of 
the favored mode are found to be negligible. 

In summary, we report a result of the measurement 
of the ratio TZdk*°i using a 711 fb _1 data sample col- 
lected by the Belle detector. We obtain lZ DK *o — 
(4.1+goi?'!) x 10~ 2 , which can be used to extract 
by combining with other observables related to the same 
dynamical parameters r$, 5s and k. Since the value 
of TZdk' is not significant, we set an upper limit of 
Rojj-.o < 0.16 (95% C.L.); this is the most stringent limit 
to date. Possible reasons for the small rs are destructive 
interference between the two D decays, destructive in- 
terference between DK*° and other DK + ir~ states, or a 
small ratio of magnitudes of amplitudes for B° — > D°K*° 
over B° -> D°K*°. 
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